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q^I Abstract. Two mechanisms by which the Quintessence scalar could enhance the relic abundance 

. of dark matter particles are discussed. This effect can have an impact on supersymmetric candidates 

' for dark matter. 
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According to the standard paradigm, a particle species goes through two main regimes 
\ during the cosmological evolution. At early times it stays in thermal equilibrium, until 

Q\ ; the particle interaction rate T remains larger than the expansion rate H. Later on, the 

particles will be so diluted by the expansion of the universe that they will not interact 
anymore and H will overcome T. The epoch at which T = H is called 'freeze-out', and 
after that time the number of particles per comoving volume for any given species will 
remain constant. This is how cold dark matter particle relics (neutralinos, for example) 
q ! are generated. 

As it can be easily understood, this scenario strongly depends on the evolution equa- 
tion for H in the early universe, which is usually assumed to be radiation-dominated. 
However, as it was already noticed some time ago [1], there is little or no evidence that 

.£h ! before Big Bang Nuclesynthesis (BBN) it was necessarily so. Non-standard scenarios 

are then worth exploring. In particular, if we imagine that for some time in the past the 
Hubble parameter was larger than usually thought (for example, due to the presence of 
some other component, in addition to radiation), then the decoupling of particle species 
would be anticipated, resulting in a net enhancement of their relic abundance. 

A natural candidate for doing that is the Quintessence scalar, which is thought to 
consitute the Dark Energy fluid dominating the present universe. In most Quintessence 
models, the cosmological scalar is assumed to become the dominant component of the 
universe after a long period of sub-domination [2], playing little or no role in the earliest 
epochs. However, this has not always to be the case, as we will show in the following. We 
will focus on the possibility of modifying the past evolution of the Hubble paramenter 
H, with the double aim of respecting all the post-BBN bounds for the expansion rate 
and of producing a mesurable enhancement of the dark matter particles relic abundance. 
In particular, we will report about two possible mechanism by which the past dynamics 
of the Quintessence scalar could significantly modify the standard evolution of the pre- 
BBN universe: an early "kination" phase and a "scalar-tensor" model. 



KINATION ENHANCEMENT 



If we imagine to add a significant fraction of scalar energy density to the background 
radiation at some time in the cosmological history, this would produce a variation in H 2 , 
depending on the scalar equation of state w*. If > w r = 1/3, the scalar energy density 
would decay more rapidly than radiation, but temporarily increase the global expansion 
rate. This possibility was explicitly considered in Ref. [3], where it was calculated that a 
huge enhancement of the relic abundance of neutralinos could be produced in this way. 
In a flat universe, a scalar field with potential V(<j>) obeys the equations 

<j)+3H<j)+dV/d<l) = ; H 2 = (d/a) 2 = Snp/3M 2 p . (1) 

For any given time during the cosmological evolution, the relative importance of the 
scalar energy density w.r.t. to matter and radiation in the total energy density p = 
pm + pr + pip depends on the initial conditions, and is constrained by the available 
cosmological data on the expansion rate and large scale structure. If the potential V(0) 
is of the runaway type, the initial stage of the scalar evolution is typically characterized 
by a period of so-called 'kination' [2] during which the scalar energy density p^, = 
2 /2 + V(0) is dominated by the kinetic contribution E k = <j> 2 /2^>V(<j)), giving = 1. 
After this initial phase, the field comes to a stop and remains nearly constant for some 
time ('freezing' phase), until it eventually reaches an attractor solution [2]. 

Then, if we modify the standard picture according to which only radiation plays a role 
in the post-inflationary era and suppose that at some time i the scalar contribution was 
small but non negligible w.r.t. radiation, then at that time the expansion rate H(i) should 
be correspondingly modified. During the kination phase the scalar to radiation energy 
density ratio evolves like p^/p r ~ aT 3(w< ''~ HV ) = aT 2 , and so the scalar contribution 
would rapidly fall off and leave room to radiation. In this way, we can respect the 
BBN bounds and at the same time keep a significant scalar contribution to the total 
energy density just few red-shifts before. The increase in the expansion rate H due to 
the additional scalar contribution would anticipate the decoupling of particle species and 
result in a net increase of the corresponding relic densities. As shown in [3], a scalar to 
radiation energy density ratio p§j p r — 0.01 at BBN would give an enhancement of the 
neutralino codensity of roughly three orders of magnitude. 

The enhancement of the relic density of neutralinos requires that at some early time 
the scalar energy density was dominating the Universe. This fact raises a problem if 
we want to identify the scalar contribution responsible for this phenomenon with the 
Quintessence field [4]. Indeed, the initial conditions must be such that the scalar energy 
density is sub-dominant at the beginning, if we want the Quintessence field to reach the 
cosmological attractor in time to be responsible for the presently observed acceleration 
of the expansion [2]. For initial conditions p<j, ^ p r we obtain instead an 'overshooting' 
behavior: the scalar field rapidly rolls down the potential and after the kination stage 
remains frozen at an energy density much smaller than the critical one. However, as 
shown in [5], more complicated dynamics are possible if we relax the hypothesis of 
considering a single uncoupled scalar. The presence of several scalars and/or of a small 
coupling with the dark matter fields could modify the dynamics in such a way that the 
attractor is reached in time even if we started in the overshooting region. 



Consider a potential of the form V(0i, §2) = M n+4 (0102) ~ n ^ 2 , with M a constant of 
dimension mass. In this case, the two fields' dynamics enlarges the range of possible 
initial conditions for obtaining a quintessential behavior today. This is due to the fact 
that the presence of more scalars allows to play with the initial conditions in the fields' 
values, while maintaining the total initial scalar energy density fixed. Doing so, it is 
possible to obtain a situation in which for a fixed p 1 ^ in the overshooting region, if we 
keep initially 0i = 02 we actually produce an overshooting behavior, while if we choose 
to start with 0i 7^ 02 ( an d the same pi n ) it is possible to reach the attractor in time. 

Suppose, instead, that the Quintessence scalar is not completely decoupled from the 
rest of the Universe. Among the possible interactions, two interesting cases are the 
following: 

V b = bH 2 f or V c = cp m ^> (2) 

If we add Vb or V c to V = M n+ (j)~ n , the potential will acquire a (time-dependent) 
minimum and the scalar field will be prevented from running freely to infinity. In 
this way, the long freezing phase that characterizes the evolution of a scalar field with 
initial conditions in the overshooting region can be avoided. A more detailed discussion, 
together with numerical examples, can be found in Ref. [4]. 



SCALAR-TENSOR ENHANCEMENT 

A different possibility arises if we consider Quintessence models in the framework of 
scalar-tensor (ST) theories of gravity (see [7] and references therein). These theories 
represent a natural framework in which massless scalars may appear in the gravitational 
sector of the theory without being phenomenologically dangerous, since they assume a 
metric coupling of matter with the scalar field, thus ensuring the equivalence principle 
and the constancy of all non-gravitational coupling constants [6]. Moreover a large class 
of these models exhibit an attractor mechanism towards GR [8], that is, the expansion 
of the Universe during the matter dominated era tends to drive the scalar fields toward a 
state where the theory becomes indistinguishable from GR. 

ST theories of gravity are defined, in the so-called 'Jordan' frame, by the action 

5g = I(b/ dAx ^[^ 2 R +4(0($)^ v ^^ v *-4t/(4>)] . (3) 

The matter fields ^ m are coupled only to the metric tensor g^ v and not to 4>, i.e. S m = 
Sm^nngfiv]- Each ST model is identified by the two functions ft)(4>) and V(<&). The 
matter energy-momentum tensor is conserved, masses and non-gravitational couplings 
are time independent, and in a locally inertial frame non gravitational physics laws take 
their usual form. Thus, the 'Jordan' frame variables g^ v and $ are also denoted as the 
'physical' ones in the literature. By means of a conformal transformation, 

g^=A 2 ((p)g^ , c£>2 = 87rM 2 A -2 ((p) (4) 



with 

a 2 ((p) = dlogA((p)/d(p = l/(4©(4>) + 6) , 



(5) 



it is possible to go the 'Einstein' frame in which the gravitational action takes the 
standard form, while matter couples to <p only through a purely metric coupling, 

S„ = S m [¥ m ,A 2 (<p)g llv ] ■ (6) 

In this frame masses and non-gravitational coupling constants are field-dependent, and 
the energy-momentum tensor of matter fields is not conserved separately, but only when 
summed with the scalar field one. On the other hand, the Einstein frame Planck mass 
is time-independent and the field equations have the simple form 

Rp V -fa v R=T* v /M* + T (lv /Ml , M 2 d 2 ( P + dV/d( P = -a(( P )T/V2. (7) 

When a (<p) = the scalar field is decoupled from ordinary matter and the ST theory 
is indistinguishable from ordinary GR. The effect of the early presence of a scalar 
field on the physical processes will come through the Jordan-frame Hubble parameter 
H = dloga/dt: 

H = H(l + a(( P )(p / )/A(( P ), (8) 

where H = dloga/dT is the Einstein frame Hubble parameter. A very attractive class 
of models is that in which the function a (<p) has a zero with a positive slope, since this 
point, corresponding to GR, is an attractive fixed point for the field equation of motion 
[8]. It was emphasized in Ref. [9] that the fixed point starts to be effective around matter- 
radiation equivalence, and that it governs the field evolution until recent epochs, when 
the Quintessence potential becomes dominant. If the latter has a run-away behavior, the 
same should be true for <x((p), so that the late-time behavior converges to GR. For this 
reason, we will consider the following choice , 

A{q>) = \+Be-P (p , a(<p) = -fiBe^ /{\ + Be~^) , (9) 

which has a run-away behavior with positive slope. 

In Ref. [7] it was calculated the effect of ST on the Jordan-frame Hubble parameter 
H at the time of WIMP decoupling, imposing on the parameters B and /3 the constraints 
coming from GR test, CMB observations and BBN. Computing the ratio H/Hgr at 
the decoupling time of a typical WIMP of mass m = 200 GeV, it was found that it is 
possible to produce an enhancement of the expansion rate up to O(10 5 ). As a further 
step, it was performed the calculation of the relic abundance of a DM WIMP with mass 
m and annihilation cross-section (<7annv)- The effect of the modified ST gravity enters 
the computation of particle physics processes (like the WIMP relic abundance) through 
the "physical" expansion rate H defined in Eq. (8). We have therefore implemented the 
standard Boltzmann equation with the modified physical Hubble parameter H: 

dY/dx=-s{o mn v){Y 2 -Y 2 q )/Hx (10) 

where x = m/T, s = (2k 2 /45) K(T) T 3 is the entropy density and Y = n/s is the WIMP 
density per comoving volume. 

A numerical solution of the Boltzmann equation Eq. (10) is shown in Fig. 1 for a 
toy-model of a DM WIMP of mass m = 50 GeV and constant annihilation cross-section 
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x = m /T 
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FIGURE 1. Numerical solution of the Boltzmann equation Eq. (10) in a ST cosmology for a toy-model 
of a DM WIMP of mass hi — 50 GeV and constant annihilation cross-section (<7 ann v) — 1 x 10 ^ GeV ^. 
The temperature evolution of the WIMP abundance Y(x) clearly shows that freeze-out is anticipated, 
since the expansion rate of the Universe is largely enhanced by the presence of the scalar field (p. At a 
value x — m/Ty a re-annihilation phase occurs and Y (x) drops to the present day value. 

(Oannv) = 1 x 10~ 7 GeV 2 . The temperature evolution of the WIMP abundance Y(x) 
clearly shows that freeze-out is anticipated, since the expansion rate of the Universe is 
largely enhanced by the presence of the scalar field (p. This effect is expected. However, 
we note that a peculiar effect emerges: when the ST theory approached GR (a fact which 
is parametrized by A (<p) — > 1 at a temperature Tm, which in our model is 0.1 GeV), H 
rapidly drops below the interaction rate T establishing a short period during which the 
already frozen WIMPs are still abundant enough to start a sizeable re-annihilation. This 
post-freeze-out "re-annihilation phase" has the effect of reducing the WIMP abundance, 
which nevertheless remains much larger than in the standard case (for further discussion 
on this aspect see [7]). The amount of increase in the relic abundance which is present 
in ST cosmology is shown in Fig. 2. The solid curve refers to an annihilation cross 
section constant in temperature, i.e. (a an nv) = a, while the dashed line stands for an 
annihilation cross section which evolves with temperature as: (c7 a nnv) = b/x. In the case 
of 5-wave annihilation the increase in relic abundance ranges from a factor of 10 up to a 
factor of 400. For a pure b/x dependence, the enhancement can be as large as 3 orders of 
magnitude. Needless to say, such potentially (very) large deviations entail new prospects 
on the WIMP characterization both for the choice of the CDM candidates and for their 
direct and indirect detection probes [7, 10]. 




FIGURE 2. Increase in the WIMP relic abundance in ST cosmology with respect to the GR case. 
The solid curve refers to an annihilation cross section constant in temperature, i.e. (tJannV) = a = 10~ 7 
GeV -2 , while the dashed line stands for an annihilation cross section which evolves with temperature as 

(cTannv) = b/x = 1(T 7 GeV- 2 /*- 
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